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ABSTRACT 

The current census of observed baryons in the local Universe is still missing a 
significant fraction of them according to standard Big-Bang nucleosynthesis. Numerical 
simulations predict that most of the missing baryons are in a hot intergalactic medium, 
which is difficult to observe through its X-ray emission or Sunyaev-Zel'dovich effect. 
We show that the next generation of X-ray satellites will be able to detect this gas 
through the X-ray absorption lines imprinted by its highly-ionized metals on the 
spectrum of a background quasar. For the metallicity typically found in intracluster 
' gas, up to 70% of the baryons produce O VIII absorption lines with an equivalent 

width ;> 0.1 eV. The spectrum of any high redshift quasar is expected to show several 
such lines per unit redshift due to intervening gaseous halos of galaxy groups. These 
lines will be detectable at a signal-to-noise ratio ;> 5 after a day of integration with 
! the future Constellation-X telescope for any of the several hundred brightest quasars 

across the sky. 
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1. Introduction 

The best estimate for the sum of the baryonic mass currently observed in stars, neutral 
hydrogen, and X-ray emitting gas in clusters of galaxies (Fukugita, Hogan, & Peebles 1998) 
falls short of the most-likely value predicted by standard Big-Bang nucleosynthesis (Buries & 
Tytler 1998). Numerical simulations suggest that most of the missing baryons reside in a hot 
intergalactic medium, which is too rarefied to be detectable through its X-ray emission or the 
Sunyaev-Zel'dovich effect (e.g., Ostriker & Cen 1996; Theuns et al. 1998). Much of this gas resides 
in the outskirts of galaxy groups or clusters, and hence is likely to possess a similar metallicity to 
that observed in the X-ray emission spectrum from the cores of these systems (Mushotzky et al. 
1996). In this paper we explore the possibility of detecting this hot intergalactic medium through 
the X-ray absorption lines imprinted by its metals on the spectrum of a background quasar. 

Detection of X-ray absorption lines from galaxy clusters had been suggested in the past as a 
means of measuring distances (Krolik & Raymond 1988; Sarazin 1989). However, because of the 
limited sensitivity of X-ray instruments at that time, these early papers restricted their attention 
to the strong lines produced in the cores of rich clusters and hence to the unlikely superposition of a 
quasar behind such cores. The next generation of X-ray telescopes, manifested at their best by the 
successor to AXAF, Constellation-X (see http:// constellation.gsfc.nasa.gov), will achieve a much 
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greater sensitivity than previously imagined. It is therefore the purpose of this paper to quantify 
the full probability distribution of equivalent widths for the intergalactic X-ray absorption lines 
imprinted on the spectrum of any random quasar. Since much of the baryonic mass is currently 
assembled into virialized objects and since the present nonlinear mass-scale is that of galaxy 
groups, most of the absorption lines would be produced by gas in the outskirts of these groups. 
In §2 we combine the Press-Schechter mass function and the singular-isothermal-sphere model 
to predict the density and temperature distribution of this gas in a standard CDM cosmology 
with cjg^-i = 0.5. We used a Bardeen et al. (1986) power spectrum with a tilted primordial 
index of n = 0.7. Because the power-spectrum is normalized on the cluster mass scale, our results 
depend only weakly on the choice of the power spectrum shape (and agree by construction with 
X-ray data for the temperature function of groups and clusters; see, e.g. Henry et al. 1995). We 
assume an Q = 1, = 0.1, h = 0.6 cosmology throughout the discussion. In §3 we compute 
the corresponding statistics of the absorption lines produced by highly-ionized metals in this gas. 
Finally, §4 summarizes the main conclusions of this work. 



2. Statistical Properties of the Gas in the Outskirts of Groups and Clusters of 

Galaxies 

Observations of the X-ray surface brightness of galaxy clusters are commonly fitted by the 
isothermal /3-model. This model assumes that the gas temperature is constant and the radial 
profile of its mass density is given by 



p gas (r,T) =po(T) 



1 + 



r ,2 



r, 



-3/9/2 



(1) 



with r c being the core radius. Typical values for the parameters of this functional fit are clustered 
around r c ~ 0.25 Mpc and (3 ~ 2/3 (Jones & Forman 1984); for specificity, we adopt these values 
throughout our discussion. The central density can be matched at large radii far from the core 
(r 3> r c ), to that of a singular isothermal profile, p gas ~ fgaskT '/ \2ir fj,m p Gr 2 ) , where / gas = f2&/fi is 
the global mass fraction of the gas (assumed to equal the universal baryon-to-total matter density 
ratio), T is the gas temperature, m p is the proton mass, and [i « 1.22 is the mean atomic weight. 
Hence we set, 

We implicitly assume that the dark matter provides the gravitational potential needed to support 
the above gas distribution, <p = —(kT/[im p ) lnp gas + const. 

At any given redshift, the comoving number density of galaxy groups or clusters with a 
temperature between T and T + dT, namely [dn(T, z)/dT]dT, can be derived by combining the 
Press-Schechter (1974) mass function, dn(M,z)/dM, with the total mass-temperature relation 
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calibrated from numerical simulations (e.g., Pen 1998), 

/ t \ 3/2 

M{T - z) = M AmTTz)) ■ (3) 

with kT$ = 4.9 keV (for $7 = 1) and M$ = 47r / o cr i t (8/i~ 1 Mpc) 3 /3 is the mass inside a sphere 
of radius 8ft. -1 Mpc filled with the cosmological critical density p cr i t . Each cluster is assigned 
a cut-off radius, -R C ut-off) which defines its total mass. Typically i? C ut-off ^ r c and so 
M(T,z) « 2kTR cut ^ oS /(fim p G). 

The column density of gas particles along the line of sight to a background source which is 
projected at an impact parameter b relative to the cluster center is given by 



iV gas (6, T) = N (T) [1 + (b/r c ) 2 \ "' " , (4) 

with Nq(T) = irpo(T)r c / /j,m p . We denote by Z(X) = Nx/Nu the fractional abundance of a heavy 
element X relative to hydrogen, and by T = N l /Nx the fraction of this element which is ionized 
i times. In collisional equilibrium, the fractional abundance of the ion X % is only a function of 
temperature (Sarazin & Bahcall 1977), i.e. T = T(T). We ignore photoionization by the X-ray 
background since we focus on the high density environments of virialized objects (in which the 
highest equivalent-width lines are expected to occur) where collisional ionization dominates. 
The collisional ionization/recombination time is typically much shorter than the Hubble time 
at the temperatures (^ 0.5 keV) and overdensities (^ 10) of interest here. Assuming a uniform 
metallicity, the column density of the ion X 1 is then 

N'ib, T) = 0.92 iV gas (6, T) Z(X) T(T) , (5) 

where 0.92 is the fraction of hydrogen atoms by number. 

The equivalent width of the absorption line produced by an ion column-density N l is given by 

W i = E^g(r), (6) 



1/2 



where 



/oo ^ 
dy [1 - exp (-T(f>(y))} , </>(y) = —= exp(-y 2 ) 



and 

2ir 2 a\ 2 c fi 



t = N l 



(E/m e c 2 )(Av x /c) ' 

Here E is the transition energy, /j is the oscillator strength for absorption, m e is the electron 
rest mass, A c = 3.86 x 10~ n cm is the electron Compton wavelength, and Avx is th e velocity 
dispersion of the element X. For a thermally broadened line Avx = \jlkT j A{X)m v (Rybicki & 
Lightman 1979), where A(X) is the atomic number of the element X. 
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3. Statistics of X-Ray Absorption Lines 



Next we would like to find the mass fraction of the hot intergalactic gas that can be probed 
by an X-ray telescope of a given sensitivity. We define W m i n to be the minimum equivalent width 
of an absorption line that can be resolved by the instrument. For a cluster temperature T, the 
column density ./V* of a given ion that would yield absorption with an equivalent width W m in 
is found by inverting equation @. The impact parameter corresponding to this density can be 
found from equations (|) and (|), namely b(W min ,T) = r c { [N^/N^W^, T)] 2 - l} 1 / 2 , where 
Nq(T) = N l (b = 0, T). Finally, the total gas mass enclosed within the cylinder defined by this 
impact parameter is given by 



M gas (> W min ,T) = 2TT 2 p (T)r 3 c 



1 + 



f b(W mhl ,T) 
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The fraction of the total baryonic mass that is traced at a redshift z through X-ray observations 
with an equivalent-width sensitivity W m { n is then 



^gas(> W n 



1 



Q h p c 



dT dn(T, Zj 



dT 



(> W min ,T) . 
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Figure 1 shows the baryonic mass fraction expressed in equation (||) at two redshifts, namely 
z = in panel (a) and z = 1 in panel (b). In both panels, the solid curve shows the result for the 
A = 18.97A absorption line of O VIII and the dashed curve refers to the A = 1.85A transition of Fe 
XXV. In this calculation we have assumed that the spectral lines are only thermally broadened, 
since the lines are generically not saturated and so turbulent broadening has a negligible influence 
on their equivalent width. The fractional abundance T(T) has been computed by interpolating 
the data from the tables of Arnaud & Rothenflug (1985) for O VIII and of Arnaud & Raymond 
(1992) for Fe XXV. We have assumed the metallicity of the gas to be 0.3^0 for iron and 0.5Z© 
for oxygen (see, e.g. Mushotzky et al. 1996; Gibson et al. 1997; Fukuzawa et al. 1998) with the 
solar metallicity abundances taken from Anders & Grevesse (1989). 

Figure la shows that at low redshifts most of the collapsed baryons can be probed by an 
instrument with a sensitivity ;> 0.1 eV to the O VIII line. The fraction is lower but still substantial 
at moderate reshifts, as shown in Figure lb. Overall, the baryonic mass probed in absorption 
extends up to an order of magnitude above the value probed by the conventional methods of 
detecting X-ray emission (Fukugita et al. 1997). The O VIII line traces more baryons than the Fe 
XXV line, since the oxygen ion exists in the temperature regime of 0.1-1 keV while the iron ion 
is abundant only at temperatures ;> 1 keV (Arnaud & Rothenflug 1985). The lower temperature 
regime is characteristic of galaxy groups (Henry et al. 1995; Burns et al. 1996), which contain 
most of the mass in the present-day Universe. 



In order to assess the feasibility of these observations it is useful to estimate the expected 
number of absorbers of a given equivalent width per unit redshift. The differential cross-section for 
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intercepting an ion column density between N l and N l + dN l within a given cluster of temperature 
T is given by [do/dN^N" 1 ) = 27rb{N i )\db/dN i \ b=b{m) . Using equations (|) and (§) we get 

da r(N\T) = { 2 ^ 2 A N h( T )?/{N l )\ < in (T) <N*< N*(T) 



dN l 1 , otherwise. 

where N^-^T) is the ion column density at b = -R cu t-off- Correspondingly, the differential 
cross-section for detecting a resonant absorption line with an equivalent width between W l and 
W i + dW l is given by [do/dW l ](W l ] = [do / dN^N 1 (W^dN 1 /dW l \ Nl=N , {Wt) . Equation (|) then 
yields, 



:(W\T)=r 



n -1 



dr 



(10) 

r=r(W i ) 



dW iK ' ; c [N^W^T)}^ TTam e c 2 \ 2 J,, 

for W l (N^ ain ) < W l < W 1 {Nq). Finally, by integrating over the cluster temperature distribution 
we derive the probability of observing an absorption line with an equivalent width between W l 
and W l + dW l per unit redshift at a redshift z, 



dW l dz 



dt 

dz 



i+ z y, (ii) 



where in the adopted = 1 cosmology the redshift derivative of cosmic path length is given by 
\cdt/dz\ = cHq 1 (1 + z)~ 5 / 2 , with Hq being the Hubble constant. Figure 2 shows this probability 
distribution at the redshifts z = [panel (a)] and z = 1 [panel (b)] for the previously mentioned 
absorption lines of O VIII (solid line) and Fe XXV (dashed line). We find that the spectrum of 
any random quasar at z = 2 is expected to show several O VIII lines with an equivalent width 
^ 0.1 eV due to intervening gaseous halos of galaxy groups. 

Finally, we would like to estimate the detectability of the predicted absorption signal for an 
X-ray telescope of a given area and spectral resolution. The maximum signal-to-noise ratio, S/N, 
is dictated by photon counting statistics. Given the characteristic surface brightness distribution 
of clusters, emission by the intracluster gas can be ignored at several core radii for a characteristic 
telescope aperture <J 1'. The required number of continuum photons from the quasar which the 
detector must receive within its spectral resolution bin of width AE around the line is (Sarazin 
1989), 



(# of photons) J ^— ) . (12) 

Correspondingly, the flux from the quasar at the line energy E needs to be 

Fe z -aqwSJ ' (13> 

where A is the collecting aperture of the telescope, Q is the quantum efficiency of the detector, 
and t is the integration time of the observation. The future Constellation-X telescope with 
its planned effective area of QA ~ 15, 000 cm 2 and resolution of AE = 1 eV could therefore 
detect an O VIII line with E = 0.654/(1 + z) keV and an equivalent-width W^in = 0.1 eV 
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at a signal-to-noise ratio S/N = 5 after t = 10 5 seconds of integration on a quasar with 
Fe ^ 2 x 10~ 12 erg cm -2 s _1 keV -1 . Based on the number counts of quasars brighter than a 
flux S in the 0.5-3.5 keV band, N(> S) ps 20 x (S/W^ 11 erg cm" 2 s~ 1 )" L5 sr" 1 (Maccacaro et al. 
1982) and the intrinsic soft X-ray spectrum of quasars Fe oc E 6 (Laor et al. 1997), we estimate 
that there should be several hundred quasars above the necessary flux threshold across the sky. 
Figure 1 therefore implies that Constellation-X will be able to probe up to 70% of the hot gas in 
the outskirts of groups or clusters of galaxies through the O VIII resonant absorption line. 

4. Conclusions 

We have investigated the possibility of using the X-ray absorption spectrum of a background 
quasar to probe the presence of heavy elements in the outer regions of groups and clusters of 
galaxies. Figure 1 demonstrates that observations of the O VIII line with an equivalent width 
sensitivity ^ 0.1 eV should be able to probe most of the baryonic material out to the virial radius 
of these systems, well beyond the regions where X-ray emission is currently detected. We estimate 
that there should be several such O VIII absorption lines per unit redshift along the line of sight 
to any quasar (cf. Fig. 2). The future Constellation-X telescope will be able to detect these lines 
with a signal-to-noise ratio ^ 5 in a day of integration for any of the few hundred brightest quasars 
accross the sky. 

Our model assumed that the gas in groups or clusters of galaxies is isothermal. The modest 
radial decline which was observed recently with ASCA for the temperature of the gas in the outer 
parts of clusters (Markevitch et al. 1998), might make absorption features by ionic species such 
as O VII more prominent than O VIII in the X-ray absorption spectra of quasars. On the other 
hand, energy input by the same supernovae that enrich the gas with metals might increase the 
temperature of X-ray halos beyond the value adopted in equation (0) (which was calibrated based 
on numerical simulations of gravitational dynamics only). 

In our calculations we have assumed a constant metallicity of the gas; a reduced metal 
abundance due to inefficient mixing in the outer parts of groups and clusters would obviously 
lower our predicted signal. If absorption will be detected for some quasars, it would be interesting 
to correlate the absorption line redshift and sky position with X-ray or optical catalogs of galaxy 
clusters. If a group or a cluster are identified at the absorption redshift, it would be possible to 
infer the impact paramter b, and the metallicity Z, using equations (|5|) and (||). Observations of 
the O VIII line around X-ray clusters at different redshifts could trace the distribution of metals in 
the intergalactic medium and possibly shed light on the epoch and nature of the metal enrichment 
process. 

We thank Richard Mushotzky and George Rybicki for useful discussions. This work was 
supported in part by the NASA grant NAG 5-7039 and by a fellowship from the university of 
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Fig. 1. — The fraction of baryons in groups or clusters of galaxies that can be probed by an X-ray 
telescope with an equivalent-width sensitivity W m i n to absorption lines. The solid curve refers to 
absorption by the strongest resonant transition of O VIII [at 0.654/(1 + z) keV, where z is the 
absorption redshift], while the dashed curve corresponds to the strongest resonant absorption by 
Fe XXV [at 6.7/(1 + *) keV]. 
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Fig. 2. — The probability for observing an absorber equivalent width W l per unit redshift per 
logarithmic equivalent-width interval as a function of the equivalent width. 



